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SUMMARY OF COMPLETED WORK 

This grant was awarded by NASA to The University of Alabama in Huntsville (UAH) 
to construct state-of-the-art, theoretical, two-component, chromospheric models for single 
stars of different spectral types and different evolutionary status. In our proposal, we 
suggested to use these models to predict the level of the “basal flux” , the observed range 
of variation of chromospheric activity for a given spectral type, and the decrease of this 
activity with stellar age. In addition, for red giants and supergiants, we also proposed 
to construct self-consistent, purely theoretical wind models, and used these models to 
investigate the origin of “dividing lines” in the H-R diagram. In the following, we describe 
our completed work. 

We have accomplished the first main goal of our proposal by constructing first purely 
theoretical, time-dependent and two-component models of stellar chromospheres.| 
The models require specifying only three basic stellar parameters, namely, the effective tem- 
perature, gravity and rotation rate, and they take into account non-magnetic and magnetic 
regions in stellar chromospheres. The non-magnetic regions are heated by acoustic waves 
generated by the turbulent convection in the stellar subphotospheric layers. The magnetic 
regions are identified with magnetic flux tubes uniformly distributed over the entire stellar 
surface and they are heated by longitudinal tube waves generated by turbulent motions in 
the subphotospheric and photospheric layers. The coverage of stellar surface by magnetic 
regions (the so-called filling factor) is estimated for a given rotation rate from an obser- 
vational relationship. The constructed models are time-dependent and are based on the 
energy balance between the amount of mechanical energy supplied by waves and radiative 
losses in strong Ca II and Mg II emission lines. To calculate the amount of wave energy 
in the non-magnetic regions, we have used the Lighthill- Stein theory for sound generation 
modified by Musielak, Rosner, Stein & Ulmschneider (1994). The acoustic wave energy 
fluxes for chromospherically active stars located in different regions of the H-R diagram 
have been computed (Ulmschneider, Theurer & Musielak 1996; Ulmschneider, Theurer, 
Musielak & Kurucz 1998). The wave energy fluxes carried by longitudinal waves along 
magnetic flux tubes have also been calculated by using both analytical and n um erical 
methods. Our analytical approach is based a theory developed by Musielak, Rosner, & 
Ulmschneider (1989) and Musielak, Rosner, Gail & Ulmschneider (1995) which allows com- 
puting the wave energy fluxes for linear waves. A numerical approach has been developed 
by Ulmschneider &: Musielak (1998) to compute the energy fluxes for nonlinear waves. Both 
methods have been used to calculate the wave energy fluxes for stars located in different re- 
gions of the HR diagram (Musielak, Rosner & Ulmschneider 1998; Ulmschneider, Musielak 
&: Fawzy 1998). Having obtained the wave energy fluxes carried by acoustic and longitu- 
dinal tube waves, we have investigated the behavior of these waves in the solar and stellar 
atmospheres (Theurer, Ulmschneider & Cuntz 1996; Theurer, Ulmschneider & Kalkofen 
1997; Sutmann, Musielak &: Ulmschneider 1998; Fawzy, Ulmschneider & Cuntz 1998) and 
use the results to construct theoretical models of stellar chromospheres (Cuntz, Ulmschnei- 
der & Musielak 1998; Buchholz, Ulmschneider & Cuntz 1998). We have used these models 
to predict theoretically the “basal flux” in late-type dwarfs and giants, and the decrease of 
the Ca II emission with decreasing stellar rotation rate. The first purely theoretical, 
two-component and time-dependent models of stellar chromospheres have been 
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constructed by Cuntz, Rammacher, Ulmschneider, Musielak & Saax (1998). 

The first theoretical, two-component chromospheric models have been constructed for 
K2 V stars with T eff = 4900 K, log g = 4.5 and the radius R = 0.8 R©. The magnetic field 
strength B 0 inside magnetic flux tubes at the photospheric level rsooo = 1 is calculated 
by assuming that the gas pressure inside and outside the tube (pi and p e , respectively) 
is given by pi/p e = 1/5, which holds for the Sun (Solanki 1996) and is extended here to 
the considered stars. This gives B 0 = 2100 G which is independent of stellar rotation. 
We assume that the magnetic flux tubes with B 0 are uniformly distributed over the entire 
stellar surface. To construct our two-component chromospheric models, we have to know 
the area of the stellar surface that is covered by magnetic and non-magnetic regions; 
this can be done by specifying the photospheric magnetic filling factor f Q . Since we are 
interested in the effects of stellar rotation rate P rot on the level of chromospheric activity, we 
have to relate B 0 f 0 to P TO t . The relationship of our choice is based on the results of Riiedi 
et al. (1997) who presented a very recent analysis of 13 high-quality optical spectra for late- 
type dwarfs (G1-K5) with different levels of activity. The derived relationship has been 
used to constrain the shape of the flux tubes (see Cuntz, Ulmschneider & Musielak 1998 for 
details). Fast rotating stars are expected to have higher coverage of flux tubes than slow 
rotating stars implying increased photospheric and chromospheric magnetic filling factors. 
As a consequence the tube spreading in fast rotating stars should be smaller. 

Having specified the structure of magnetic regions, the wave calculations are performed 
by solving a set of time-dependent, nonlinear and ideal MHD equations for thin magnetic 
flux tubes. The equations are solved using the method of characteristics (Fawzy, Ulm- 
schneider & Cuntz 1998) and the wave energy fluxes obtained by Ulmschneider, Musielak 
and Fawzy (1998) are used. In non-magnetic regions, the wave calculations Eire performed 
by solving time-dependent and nonlinear hydrodynamic equations for acoustic waves and 
the wave energy fluxes computed by Ulmschneider, Theurer & Musielak 1996) are used. 
In both cases of magnetic and non-magnetic regions, only monochromatic waves are con- 
sidered. In addition, the models for both regions are computed separately, which means 
that there is no mutual interaction between both components and there is no energy leak- 
age from one component to the other. Simultaneously with these wave calculations, we 
perform the computations of radiative losses in H“ continuum and in Mg II and Ca II 
emission lines by solving the appropriate radiative transfer equations together with the 
statistical equilibrium equations for NLTE populations. The time-dependent energy bal- 
ance between the wave dissipation due to the formation of shocks and the emitted radiative 
losses is calculated at different heights in the atmosphere and it is used to determine the 
local values of temperature, density and pressure. The time-averaged chromospheric mod- 
els of magnetic and non-magnetic regions are calculated by averaging over a timespan of 
many wave periods. Since the heating in magnetic regions is much more efficient than in 
non-magnetic regions, the increase in the level of chromospheric activity in stars of the 
same spectral type is obtained by increasing the coverage of the stellar surface by magnetic 
regions. To achieve this, higher magnetic filling factors axe required, which imply higher 
stellax rotation rates — a trend fully consistent with the observational data. 

The first theoretical, two-component chromospheric models have been constructed for 
K2 V stars with T e ff = 4900 K, log g = 4.5 and the radius R = 0.8 R©. The magnetic field 
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strength B 0 inside magnetic flux tubes at the photospheric level r 50 oo = 1 is calculated 
by assuming that the gas pressure inside and outside the tube (pi and p e , respectively) 
is given by Pi / Pe = 1/5, which holds for the Sun (Solanki 1996) and is extended here to 
the considered stars. This gives B 0 = 2100 G which is independent of stellar rotation. 
We assume that the magnetic flux tubes with B 0 are uniformly distributed over the entire 
stellar surface. To construct our two-component chromospheric models, we have to know 
the area of the stellar surface that is covered by magnetic and non- magnetic regions; 
this can be done by specifying the photospheric magnetic filling factor f 0 . Since we are 
interested in the effects of stellar rotation rate P ro t on the level of chromospheric activity, we 
have to relate B 0 f 0 to P TO t . The relationship of our choice is based on the results of Riiedi 
et al. (1997) who presented a very recent analysis of 13 high-quality optical spectra for late- 
type dwarfs (Gl— K5) with different levels of activity. The derived relationship has been 
used to constrain the shape of the flux tubes (see Cuntz, Ulmschneider & Musielak 1998 for 
details). Fast rotating stars are expected to have higher coverage of flux tubes than slow 
rotating stars implying increased photospheric and chromospheric magnetic filling factors. 
As a consequence the tube spreading in fast rotating stars should be smaller. 


Having specified the structure of magnetic regions, the wave calculations are performed 
by solving a set of time-dependent, nonlinear and ideal MHD equations for thin magnetic 
flux tubes. The equations are solved using the method of characteristics (see Herbold et 
al. 1985; Rammacher & Ulmschneider 1989; Fawzy, Ulmschneider & Cuntz 1998) and 
the wave energy fluxes obtained by Ulmschneider, Musielak and Fawzy (1998) are used. 
In non-magnetic regions, the wave calculations are performed by solving time-dependent 
and nonlinear hydrodynamic equations for acoustic waves (see Ulmschneider, Muchmore & 
Kalkofen 1987; Rammacher & Ulmschneider 1992) and the wave energy fluxes computed 
by Ulmschneider, Theurer & Musielak 1996) are used. In both cases of magnetic and 
non-magnetic regions, only monochromatic waves are considered. In addition, the models 
for both regions axe computed separately, which means that there is no mutual interaction 
between both components and there is no energy leakage from one component to the other. 
Simultanously with these wave calculations, we perform the computations of radiative 
losses in H continuum and in Mg II and Ca II emission lines by solving the appropriate 
radiative transfer equations together with the statistical equilibrium equations for NLTE 
populations (Ulmschneider, Muchmore & Kalkofen 1987; Rammacher & Cuntz 1991). The 
time- dependent energy balance between the wave dissipation due to the formation of shocks 
and the emitted radiative losses is calculated at different heights in the atmosphere and it is 
used to determine the local values of temperature, density and pressure. The time-averaged 
chromospheric models of magnetic and non-magnetic regions are calculated by averaging 
over a timespan of many wave periods. Since the heating in magnetic regions is much more 
efficient than in non-magnetic regions, the increase in the level of chromospheric activity in 
stars of the same spectral type is obtained by increasing the coverage of the stellar surface 
by magnetic regions. To achieve this, higher magnetic filling factors are required, which 
imply higher stellar rotation rates — a trend fully consistent with the observational data. 
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Figure 1. The decrease of Ca II fluxes with decreasing stellar rotation rate, (a) Theoretical 
Ca II fluxes (squares) calculated by using our theoretical, two-component chromospheric 
models are directly compared to the observational data (triangles) for a set of stars with 
spectral type between KO V and K3 V. (b) The comparison between the observations at 
low activity phases of the stars (triangles) and the theoretical basal flux limit is shown 
(after Cuntz, Rammacher, Ulmschneider, Musielak & Saar 1998). 

The computed purely theoretical chromospheric models are then used to calculate the 
emergent Ca II H-fK line fluxes and profiles, and compare them with the observational 
results (see Fig. 1). Note that for the purpose of calculating these emergent fluxes and 
profiles, the magnetic and non-magnetic regions must be composed together (see C un tz, 
Rammacher, Ulmschneider, Musielak & Saar 1998). The results presented in panel (a) of 
Figure 1 show the decrease of Ca II flux with decreasing stellar rotation rate. Theoretically 
calculated Ca II fluxes (denoted by squares) axe directly compare to the observational 
data (denoted by triangles) for a set of stars with spectral type between KO V and K3 V. 
In panel (b) of the same figure the comparison between low activity phases of the stars 
obtained via long-term monitoring programs (triangles) and the theoretical basal flux limit 
(Buchholz, Ulmschneider & Cuntz 1998) is shown. In both cases the agreement between 
the theoretical results and the observational data is very good, which shows that our models 
are able to reproduce empirical chromospheric emission - stellar rotation relations as well 
as chromospheric basal flux limits for stars of low activity. This gives us confidence that 
our two-component chromospheric models are consistent with observations and 
that we are on the right track to explain the origin and evolution of stellar chromospheres. 

We have also accomplished our second main goal of the proposed research by con- 
structing theoretical, time-dependent and self-consistent wind models based on 
the momentum deposition by nonlinear Alfven waves. The full set of single-fluid MHD 
equations with the background flow has been solved by using a modified version of the 
ZEUS MHD code. The constructed wind models are radially symmetric with the magnetic 
field decreasing radially and the initial outflow is described by the standard Parker wind 
solution. In contrast to previous studies, no assumptions regarding wave linearity, wave 
damping, and wave-flow interaction are made; the models thus naturally account for the 
backreaction of the wind on the waves as well as for the nonlinear interaction between 
different types of MHD waves. The models have been already used to explain the origin 
of fast speed streams in solar coronal holes (see Ong, Musielak, Rosner, Suess & Sulkanen 
1997). The obtained results clearly demonstrate that the momentum deposition by Alfven 
waves in the solar wind can be sufficient to explain the origin of fast stream components of 
the solar wind. The range of wave amplitudes required to obtain the desired result seems to 
be in good agreement with recent observations. After obtaining these encouraging results, 
we intended to use the developed code to construct self-consistent models of cool massive 
winds observed from late-type giants and supergiants. The work was supposed to be done 
by Mr. Ong as part of his Ph.D. research. Unfortunately, Mr. Ong, who had been working 
on this project for the last two years and was partially supported by this grant, decided to 
accept a position in the computer industry and discontinued working on his dissertation. 
In the following, we briefly describe these results and then show how to apply them to cool 
massive winds observed from the evolved stars. 
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Figure 2. Time-averaged wind velocity profiles for four different solar wind models. In 
each model, the initial flow is assumed to be the thermally-driven (Parker) wind. Alfven 
waves with the same period (1 hour) but different amplitudes, A, are generated at the base 
of each model (corresponding to the base of solar coronal holes) and allowed to interact with 
the initial flow for 278 hours. After this time, the resulting time-averaged wind velocity 
profiles are calculated. The obtained results show that the presence of Alfven waves with 
amplitudes ranging from 40 km s -1 to 75 km s -1 is required in the solar wind in order to 
explain the origin of fast streams (after Ong, Musielak, Rosner, Suess k Su lkan en 1997). 

To calculate our wind models, we start with a standard Parker model of thermally 
driven winds (Parker 1958) and treat it as the initial flow. The structure of an atmosphere 
computed with this initial flow is then perturbed by radially-propagating toroidal Alfven 
waves of a finite amplitude. The interaction between these waves and the flow is treated 
self-consistently, e.g., the structure of the wind and atmosphere is modified by the waves 
and, in turn, the wave behavior is influenced by this newly modified structure. The physical 
parameters describing the background medium are recorded (as they fluctuate on time 
scales of the Alfven wave period), and then time-averaged over the Alfven wave crossing 
time in order to obtain estimates for the corresponding parameter values of the flow on 
time scales long when compared to the Alfven wave period. In order to perform these 
calculations, we modified the original ZEUS code (Stone k Norman 1992a, b), and use 
it to solve the full set of ideal single fluid compressible MHD equations in one-dimension; 
we consider a single magnetic field line along which the field decreases radially. At the 
lower boundary of our computational domain the field line is perturbed in such a way that 
toroidal (d/d<f> = 0) Alfven waves are continuously generated. The upper boundary is 
transparent, which means that the wind and waves freely leave the computational domain. 
Since we are solving the single-fluid MHD equations without any further simplifications, 
processes such as wave reflection, nonlinear coupling between the various MHD waves, and 
shock formation are automatically accounted for. 

Our results are presented in Figure 2, which shows the distribution of the time-averaged 
wind velocity with distance. The following physical parameters at the base of the solar 
corona (r = 1 Rq) are used: the strength of the magnetic field B 0 = 10 G, temperature 
T 0 = 1.0 x 10 6 K and density p Q = 1.3 x 10 -15 . The obtained results clearly show that 
Alfven wave amplitudes ranging from 40 km s -1 to 75 km s -1 at the coronal base of the 
solar wind to explain the acceleration of fast speed streams in the solar wind. This range 
of required Alfven wave amplitudes is in agreement with observations of line broadening of 
Mg coronal emission lines (Hassler et al. 1990), with more recent Ulysses observations (e.g., 
Phillips et al. 1995), and with very recent independently performed numerical simulations 
by Ofman k Davila (1998). 

In addition to these major accomplishments, we have completed our analytical (Musiel ak , B 
Rosner k Ulmschneider 1998a) and numerical (Huang, Musielak k Ulmschneider 1995) in- 
vestigation of the generation of transverse magnetic tube waves in stellar convection zones; 
note that these waves have not yet been incorporated in our stellar chromosphere models 
described above. The results of these studies will be used to calculate the wave energy fluxes 
carried by these waves in stars located in different regions of the HR diagram (Musielak, 
Rosner k Ulmschneider 1998b; Musielak, Ulmschneider k Fazwy 1998). We have also 
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investigated the efficiency of the energy transfer along magnetic structures (Huang 1996; 
Wu, Xiao, Musielak Sc Suess 1996; Ziegler Sc Ulmschneider 1997a, b; Huang, Musielak Sc 
Ulmschneider 1998a,b). The primary goal of this study is to investigate the validity of the 
thin flux tube approximation used in our chromospheric models. The relevant work to our 
wind modeling project has been done by Krogulec Sc Musielak (1998) and by R. Rosner 
(in several papers published with different collaborators - see Section I.C) who have inves- 
tigated the behavior of linear and nonlinear Alfven waves in solar coronal holes and stellar 
atmospheres and the effects caused by these waves on the background medium. M. Cuntz 
has also worked on ab-initio models for outer atmospheric flows in a Ori (Cuntz 1997) 
and the interpretation of chromospheric velocity fields in various late-type (super- )giants 
considering proposed heating mechanisms (Harper, Bennett &: Cuntz 1998). 

The results described above have been obtained by the P.I. (Dr. Z. E. Musielak), Co-I’s 
(Drs. R. Rosner and P. Ulmschneider), one senior research associate (Dr. M. Cuntz, who 
joined UAH in Jan. 1996), one junior research associate (Dr. P. Huang, who joined UAH 
in Jan. 1996 and left in May 1996 to work for industry), and two graduate students in 
physics; one of them, Mr. K. K. Ong, has recently left UAH to work for industry without 
completing his Ph.D. Dr. Cuntz has devoted all his time to work on the project and has 
been fully supported by the grant. He has been closely working with the P.I. and Dr. 
Ulmschneider, and with Mr. Ong while he was at UAH. He also worked with Dr. Huang 
during her stay at UAH. The Co-I’s visited UAH several times and the P.I. spent some 
time at both the University of Chicago and the University of Heidelberg. In Fall of 1996, 
Mr. Ong spent three months at the Univeristy of Chicago working with Dr. Rosner on 
construction of self-consistent and time-dependent stellar wind models. Finally, the P.I. 
has been working on several problems directly related to the project during the regular 
academic year when his salary is fully paid by UAH and during the academic year 1997/98 
when he was on sabbatical leave. 
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